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ABSTRACT 

This  paper  explains  the  experiment  that  HF 
skywave  radar  probed  the  artificial 
ionospheric  hole,  which  is  caused  by  the 
flames  of  the  rocket  vertically  launched  and 
penetrated  the  ionosphere.  After  the  rocket 
had  passed  through  the  ionosphere,  the 
minimum  time-delay  Pmin-f  on  the  backscatter 
ionograms  obviously  appeared  with  wave  and 
focusing  stripes  resulting  from  the  irregular 
structure  etc.  The  results  indicated  that  there 
was  a  low  electron  density  zone,  the 
artificially  created  hole  in  the  ionosphere 
along  the  propagation  path.  Under  the 
asymmetry  quasi-cosine  ionospheric  hole 
model,  the  experimental  Pmin-f  was  simulated 
with  the  technique  of  ray  tracing.  It  was 
deduced  that  the  range  size  of  the  hole  in 
ionosphere  along  the  radar  beams  was  some 
573  km  and  the  critical  frequency  of  the  center 
was  2.6  MHz  lower  than  the  background 
critical  frequency  (12MHz).  The  ionospheric 
environment  350km  away  from  the  launching 
site  was  disturbed  and  the  propagation 
velocity  of  the  ionospheric  disturbance  was 
about  50m/s. 

Key  words:  HF  backscatter  propagation, 
artificial  ionospheric  hole,  HF  sky-wave  radar, 
2-  D  ray  tracing 

I  INTRODUCTION 

On  14  May  1973,  the  US  Scientists 
observed  the  ionospheric  effects  of  the 
launching  of  the  Space  Laboratory  with  the 


launching  rocket  of  Saturn  V  by  Faraday 
effect  technique  along  the  Atlantic  Coast  in 
the  northern  part  of  USA.  They  found  that 
there  were  large  quantity  of  electron 
compound  in  the  ionosphere  due  to  the  spurt 
flames  of  the  rocket  which  contained  large 
quantity  of  hydrogen  and  water  molecules 
thus  the  artificial  ionospheric  hole  was  formed 
which  had  the  electron  density  much  lower 
than  that  of  the  background  in  range  of  about 
1000km[1].  In  1980s  three  observations  with 
vertical  ionosphere  sounding  for  the  vertical 
launching  of  large  sized  rockets  penetrating 
the  ionosphere  in  China  showed  that  there 
were  obvious  descend  of  ionospheric  critical 
frequency  with  the  abnormality  of  the 
echo[2’6]. 

We  first  observed  the  ionospheric  effect 
with  the  HF  skywave  backscattering  radar 
after  the  launching  of  a  large  sized  rocket 
penetrating  the  ionosphere  and  obtained  the 
expected  results.  This  was  the  first  successful 
observation  of  the  artificial  ionospheric  hole 
with  such  radar.  Under  the  asymmetry 
quasi-cosine  ionospheric  hole  model,  the 
experimental  Pmin-f  curve  was  simulated  to  fit 
with  the  technique  of  ray  tracing.  The  range 
size  of  the  hole  in  ionosphere  and  the 
propagation  velocity  of  the  ionospheric 
disturbance  had  also  been  deduced. 

II  DESIGN  AND  BASIS  OF  THE 
EXPERIMENT 

We  have  done  the  analysis  with  the 
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technique  of  ray  tracing  for  the  affection  of  the 
minimum  time  delay  (Pmin-f)  of  the 
backscattering  ionogram  caused  by  the 
disturbed  ionosphere  [3,4]  and  found  that  the 
bend  fluctuation  of  the  Pmin-f  in  the 
ionospheric  backscattering  ionogram  could  be 
caused  by  the  small  scaled  disturbance  of  the 
electron  density  of  ionosphere  which  shape 
was  relevant  to  the  horizontal  gradient  of  the 
ionosphere  and  its  changes.  This  had  been 
verified  by  the  observation  of  the  annular  solar 
eclipse  effect  happened  on  23  Sept.  1987  with 
the  HF  Sky-wave  Radar.  So,  if  the  there  was 
large  quantity  electron  compound  emerging  in 
the  ionosphere  due  to  the  rocket  flame  thus 
causes  the  artificial  zone  of  low  electron 
density,  then  there  should  be  bend  fluctuation 
of  the  Pmin-f  in  the  ionospheric  backscattering 
ionogram.  When  the  hole  was  formed,  there 
should  be  obvious  horizontal  positive  and 
negative  ionospheric  sectional  gradient 
changes  and  the  corresponding  turning  point 
at  3Pmin/3f  =^>  0  in  ionospheric  backscattering 
ionogram.  So,  the  frequency  sweep 
backscattering  ionogram  was  taken  as  the 
basic  data  for  analysis. 

The  beam  azimuth  of  the  probing  radar  was 
offset  about  200km  northwest  from  the 
launching  site  in  order  try  to  obtain  the 
information  of  scale,  existent  period  and  the 
horizontal  disturbance  velocity  of  the  artificial 
ionospheric  hole.  So,  the  beam  was  not 
pointed  to  the  center  of  the  launching  site  but 
intersected  with  edge  of  the  imaginary  hole. 

The  HF  radar  was  working  in  the 
auto-frequency  sweep  mode  with  the 
sweeping  range  from  5-25MHz.  The 
experimental  record  was  the  backscattering 
frequency  sweeping  ionogram  P~f  but  the 
curves  of  minimum  time  delay  Pmin-f  have 
been  analysed  as  a  main  data. 

On  the  launching  day,  from  15h30mm  to 
19h00mm  the  ionosphere  was  observed  twice 
every  ten  minutes,  other  time  once  every  30 


minutes. 

Ill  THE  OBSERVED  DATA 

A.  Fluctuated  Pmin-f  of  the  HF  backscatter 
ionogram 

The  launching  time  of  rocket  was  16h10mm 
and  the  observed  backscattering  ionogram 
kept  normal  without  fluctuation  till  17h05mm. 
From  17h  05min  to  17h50min,  there  appeared  a 
clear  bending  in  the  minimum  time  delay 
Pmin-f  between  13~20MHz  frequency  and  in 
the  distance  region  of  700-  1700km.  At 
17h23mm  the  obvious  fluctuation  reached  the 
maximum.  And  it  was  recovered  to  normal  till 
18h00m  to  the  shape  of  sunset  and  night.  Fig.  1 
shows  the  variation  of  time  sequence  of  Pmin-f 
in  the  backscattering  ionogram  after  the 
launching  of  the  rocket  at  16h04mm'  the  dot 
line  showing  the  non-fluctuated  minimum 
time  delay  as  reference. 


Fig.l  Variation  of  time  sequence  of  Pmin-f  in 
backscattering  ionogram 

B.  Abnormal  additional  tracks  in  ionogram 

6mm  after  the  launching  of  the  rocket, 
namely  at  16h10mm,  except  the  quasi-vertical 
reflection  echo  in  the  quasi-vertical  reflection 
zone  in  the  ionogram,  there  emerged  some 
abnormal  additional  tracks  (it  is  also  called 
branch  tracks)  which  was  a  little  bit  lower 
than  the  quasi-reflecting  track  in  F  layer  with 
the  virtual  height  of  about  260km,  and  the 
reflecting  frequency  from  9.4—  13.2MHz  point 
out  by  the  arrow  in  Fig.  2a.  The  strange  tracks 


589 


had  appeared  three  times  in  the  ionogram 
obtained  in  the  vertical  sounding  station  50 
km  away  from  the  launching  site  when  the 
rocket  penetrating  the  ionosphere  [2].  and  the 
time  were  154s,  164s,  190s  respectively. 


Fig.  2  a)  The  abnormal  tracks  in  ionogram  at  16h10min, 
b)  the  focus  stripe  of  the  irregular  in  ionogram  at 
16h23min  (point  out  by  the  arrow),  the  white  line  is 
the  calculated  value  of  Pmin-f  with  the  ionospheric 
hole  model 

C.  Focus-strips  of  irregular  in  ionogram 

After  17h14mm  there  were  focus-strips  of  the 
irregular  in  the  backscattering  ionogram  along 
the  oblique  range  of  1000  km.  The  focus  strips 
of  the  irregular  in  the  backscattering  at 
16h23mm  as  shown  by  the  arrow  in  Fig.2b 
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Where  vi,  v2,  v3  are  the  geocentric  angle 
positions  (rad)  of  the  front,  center  and  the 
back  calculated  from  the  radar  station.  fc  (0) 
the  critical  frequency(MHz)  of  the  top  side  of 
the  radar  station,  x  is  the  background  negative 
gradient  of  the  critical  frequency  (MHz/rad). 
(k\+Xi)!2  is  the  disturbed  range  (km)  of  the 
ionosphere  in  the  reflecting  area  along  the 
direction  of  radar  beam.  Afc  is  the  fluctuation 
amplitude  of  the  critical  frequency  of  the  hole 
(MHz/rad).  rm,  rb,  are  the  height  of  maximum 
electron  density  and  the  bottom  height  in 
ionosphere  calculated  from  the  geo-center 
respectively.  The  model  was  at  about  0  =  vi, 
v2,  v3  and  fc  (0)  and  the  value  of  first  order 
differential  coefficient  respectively  as: 


IV  ANALYSIS  AND  DISCUSSION 


A.  Model  of  artificial  ionospheric  hole 

In  order  to  study  the  artificial  ionospheric 
hole  with  the  technique  of  fitting  of 
two-dimensional  ray  tracing,  a  model  of 
asymmetric  quasi-cosine  hole  of  the  variation 
of  electron  density  with  height  and  the 
distance  of  geocentric  angle  0  was  considered, 
if  it  could  be  simply  expressed  with  the 
critical  frequency,  it  should  be 
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de  de 


=  x; 


It  can  be  seen  from  the  above  that,  fc  (0) 
and  its  differential  coefficient  at  the  joint  point 
0  =  vi,  v2,  v3  are  successive,  because  the 
vertical  distribution  of  electron  density  is  also 
successive. 


where, 
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B.  Parameters  of  artificial  ionospheric 
hole 

With  the  technique  of  two-dimension  ray 
tracing  and  in  the  model  of  the  above  stated 
ionospheric  model  of  asymmetric  quasi-cosine 
hole,  the  Pmin-f  of  the  ionogram  at  17h30min  was 
simulated.  It  shows  that  model  below  with  the 
parameters  of  the  ionospheric  hole,  the 
calculated  value  of  Pmin-f  has  good  agreement 
with  that  of  the  experimental  data.  The 
parameters  of  the  ionospheric  hole  model 
were:  rb=6570km,  rm=6700km,  /c(0) 

=12.4MHz,  x—  -5 MHz/rad,  A/c  =  -1.3MHz/rad, 
A,i=764.4km  A,2=382.2km  Vi=0.035rad,  V2 
=0.095rad,  V3=0.125rad,  so,  the  maximum 
disturbing  amplitude  of  the  critical  frequency 
was  2.6MHz  lower  than  the  background 
critical  frequency  and  the  width  of  the 
disturbed  zone  along  the  detecting  path  was 
573.3km. 

Fig.3  shows  the  profile  of  the  electron 
density  of  the  ionosphere  of  the  model  F  (the 
variation  of  the  ionospheric  electron  density 
with  0)  and  the  trace  of  the  minimum  time 
delay  P^n  vs  all  frequencies  atl7h30mm.  In 
fig.3  ray  a,  b,  ...  j  is  corresponding  Pmin  at 
frequency  13,  14,  ...  24MHz respectively. 


Fig.3  At  17h30min  the  Pmin  of  model  of  asymmetric 
quasi-cosine  hole  simulate  by  ray  tracing 

Because  the  ray  was  offset  200  km  from  the 
launching  site,  the  above  parameters  are  not 
the  parameters  of  the  artificial  ionospheric 
hole,  but  that  of  the  tangential  section  of  the 
radar  ray  around  the  edge  of  the  hole. 
Therefore,  according  to  the  relative  geometric 
positions  of  the  radar  ray  and  the  launching 
site,  it  was  deducted  that  ionosphere  350  km 


away  from  the  launching  site  had  also  been 
disturbed  by  the  rocket  flame.  In  other  words 
if  assume  the  artificial  ionospheric  hole  is  a 
circular  hole,  then  the  size  of  hole  was 
deduced  about  700km.  It  is  a  pity  that  the 
parameters  of  the  ionospheric  hole  is  not 
approved  because  of  the  lack  of  the  vertical 
sounding  data  from  the  directions  of  eastward 
or  northeastward  from  the  launching  site. 

C.  Existent  duration  of  hole  and 
propagation  velocity  of  the  disturbance 

The  HF  backscattering  ionogram  shows  that 
the  ionosphere  along  the  direction  of  radar 
beam  (200  km  away  from  the  launching  site) 
was  disturbed  60  minutes  after  the  launching 
and  the  deducted  propagation  velocity  of  the 
disturbance  is  50.5m/s.  The  apparent 
ionospheric  disturbance  along  the  direction  of 
radar  beam  lasted  more  than  40  minutes. 

Considering  the  fact  that  the  ionospheric 
disturbance  above  the  launching  site  lasted 
two  hours  and  ten  minutes,  which  was  shown 
by  the  data  obtained  by  a  vertical  ionospheric 
sounder  near  the  launching  site,  it  can  be 
deducted  that  the  Maximum  time  of  the 
disturbed  zone  caused  by  the  rocket  flames 
17h23mm  namely,  79  minutes  after  the 
launching  of  the  rocket  and  the  ionosphere 
350  km  away  from  the  launching  site  was  also 
disturbed.  Later,  the  edge  of  the  Maximum 
disturbed  zone  began  to  recombine  gradually, 
and  the  ionosphere  along  the  direction  of  radar 
beam  recovered  after  27  minutes.  Then  the 
disturbed  area  continually  narrowing  till 
1 8h  1 6mm  the  recovery  of  the  ionosphere  above 
the  launching  site.  The  whole  course  lasted 
132minutes. 

D.  Explanation  of  the  energy  focus  strips 
of  irregular 

There  was  obvious  horizontal  gradient 
along  the  direction  of  radar  beam  because  of 
the  existence  of  the  low  electron  density  area. 
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At  this  time,  usually,  the  uniqueness 
characteristics  of  the  extremum  of  the  time 
delay  P  vs  angle  in  the  backscattering 
ionogram  obtained  from  the  stratified 
homogeneity  model,  had  been  damaged  to 
form  two  multi  extremum  of  the  complex 
curves  or  more  [5].  Therefore,  the  time  delay 
corresponding  to  the  points  of  extremum 
would  not  be  the  minimum  time  delay.  And 
near  the  points  of  extremum  dP/da  =  0,  the 
flux  of  the  ray  within  the  angle  or  the  unit 
time  increased,  so  the  energy  outside  of  these 
time  delays  focused.  From  this,  it  is  known 
that  among  these  extremums,  the  one 
corresponding  to  the  minimum  time  delay  was 
the  minimum  time  delay.  Otherwise,  it  was  the 
time  delay  of  focused  energy.  Time  delays  of 
the  focused  energy  at  different  frequencies 
formed  the  focused  strips  on  the  ionogram.  So, 
focused  strips  were  the  important 
characteristics  of  the  backscattered  ionogram 
of  the  disturbed  ionosphere. 

E.  The  Target  Characteristic  of  the  Rocket 

Researches  show  that  the  radar  cross 
section  of  the  rocket  body  is  determined  by  its 
physical  dimension  which  exists  in  the 
resonance  range  in  the  HF  band.  But  the  radar 
cross  section  of  the  rocket  flame  varies  with 
the  height  in  HF  band  and  the  maximum  value 
of  RCS  reaches  40dB.  The  spectrum  analysis 
shows  that  echo  of  the  rocket  has  the 
characteristics  of  double  peak  spread,  namely 
the  combination  of  spectrums  of  the  rocket 
body  and  that  of  the  flame.  The  latter  usually 
has  wider  spreading  in  spectrum.  If  the  rocket 
turns  off  without  flame,  it  is  very  difficult  for 
the  HF  Sky-wave  radar  to  find  the  rocket 
target.  However,  the  artificial  ionospheric  hole 
due  to  the  ionospheric  electron  combined  by 
the  large  quantity  of  neutral  molecular  caused 
by  the  rocket  flame  that  makes  descending  of 
the  electron  density  can  be  fond  by  the  HF 
sky-wave  radar.  Therefore  it  is  possible  to 


fond  and  trace  the  target  like  rocket  from  the 
launching  of  the  rocket  till  the  rocket  enter  the 
ionosphere  by  the  HF  sky-wave  radar. 

V  CONCLUSIONS 

The  analysis  of  the  experimental  data  shows 
that  it  is  feasible  to  probe  the  artificial 
ionospheric  hole  with  HF  Backscatter  Radar. 
As  the  first  experiment  for  such  purpose  it  was 
also  a  successful  one.  And  it  was  also  feasible 
and  reasonable  for  the  deduction  of  the 
parameters  of  the  low  electron  density  zone  of 
the  ionosphere  with  the  technique  of  fitting  of 
two-dimensional  ray  tracing  and  the 
experimental  data  under  the  model  of 
asymmetric  quasi-cosine  hole. 

The  target  characteristic  of  the  rocket  shows 
that  it  is  possible  to  fond  and  traced  the  target 
like  rocket  from  the  launching  of  the  rocket 
till  the  rocket  to  enter  the  ionosphere  by  the 
HF  sky-wave  radar. 
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